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Abstract 

^ : 

We use a model previously formulated based on the double exchange mech- 
! anism and diagonal disorder to calculate magnetization and conductivity for 

^ ! 

^ , Lai_2:Sra;Mn03 type crystals as a function of temperature. The model repre- 

00 \ sents each Mn^^ ion by a spin S=l/2, on which an electron can be added to 

00 
O 
^ . 

Q> , actions: exchange J, and Coulomb U, and, to represent in a simple way the 



produce Mn'^+. We include a hopping energy t, two strong intratomic inter- 



^ \ effects of disorder, a Lorentzian distribution of diagonal energies of width F 

s, 

, at the Mn sites. 

c : 

O . In the strong coupling limit, J, C/ >> i,r, the model results can be ex- 



pressed in terms of t and T. 

We use the results of the model to draw "phase diagrams" that separate 
ferromagnetic from paramagnetic states and also "insulating " states where 
the Fermi level falls in a region of localized states from "metallic " where the 
Fermi level falls in a region of extended states. 

Finally, assuming that particles in extended states make the largest contri- 
bution to conductivity, we calculate the resistivity for different concentrations 
and magnetic fields and compare with experiment. 
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We conclude that for the model can be used successfully to represent the 
transport properties of the systems under consideration. 



Typeset using REVT^ 



I. INTRODUCTION 



The discovery of "colossal "magnetoresistance in Lai_2;Sra;Mn03 type compounds ^ 
and its relation to possible applications to magnetoresistance (MR) devices has attracted 
the attention of the physics community in the last times. 

Before the discovery of "colossal " MR, the earlier studies by Jonker and Van Santen 
established a temperature-doping phase diagram separating metallic ferromagnetic from 
insulating antiferromagnetic phases. Zener |^ proposed a "Double Exchange " (DE) mech- 
anism to understand the phase diagram of these compounds and the intimate link between 
their magnetic and transport properties. This DE mechanism was used by Anderson and 
Hasegawa |^ to calculate the ferromagnetic interaction between two magnetic ions, and by 
de Gennes to propose canting states for the weakly doped compounds. Kubo and Ohata 
used a spin wave approach to study the temperature dependence of the resistivity at 
temperatures well below the critical temperature and a mean field approximation at T near 
Tc. . Mazzaferro, Balseiro and Alascio used a mixed valence approach similar to that 
devised for TmSe combining DE with the effect of doping to propose the possibility of a 
metal insulator transition in these compounds. 

More recently, a wealth of experimental results have been obtained on the transport, 
optical, spectroscopic and thermal properties of these materials under the effects of external 
magnetic fields and pressures 

From the theoretical point of view, Furukawa ^ has shown that DE is essential to the 
theory of these phenomena, while Millis et al. |jTo| have argued that DE alone is not sufficient 
to describe the properties of some of the alloys under consideration and have proposed that 
polaronic effects play an important role. In a previous paper we have explored a semi- 



phenomenological model that includes the effect of disorder in the transport properties [11 



Miiller-Hartmann and Hirsch [|T2| have pointed out that a new phase appears in the proper 
derivation of the effective hopping, but have not studied its effect in the physical properties 
of the systems under consideration. 
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Quantitative comparison between calculated and measured resistivities is scarce, with 



some exceptions as is the case in |T3| where the connection between magnetization and 
resistivity is clearly shown. 



In our previous paper |Tl| we treat the Hamiltonian proposed for these systems using an 
alloy analogy approximation to the exchange terms and including the effects of disorder by 
introducing a continuous distribution of the diagonal site energies. 

Here we continue that treatment by proposing a Free Energy that allows to determine the 
magnetization as a function of temperature. We then proceed to find the Fermi energy and 
the mobility edge (ME) as functions of temperature. Finally, assuming that the conductiv- 
ity is dominated by particles occupying extended states, we draw resistivity-vs-temperature 
curves. We compare our results with experiments in single crystals of Lai_^Sra,Mn03 re- 



ported by Tokura et al. [0] for four different values of the concentration of Sr and four 
different values of the magnetic field finding that the model allows a clear description of the 
experimental results. 

In Section II we describe the results of our previous paper and the approximations made 
to obtain the conductivity, and Section III is devoted to comparison with experiment and 
discussion of the results. 

II. MODEL 



In our previous paper [11] we consider a simplified model Hamiltonian given by 



<i,j>,fJ. 



+ UY^ cl^Ci^cl^Cii -JY^S.a, (1) 

i i 

where c|^, Cj^ creates and destroys an itinerant electron with spin /i at site i respectively, 
ej are the site diagonal energies that depend on the site neighborhood. For simplicity we 
consider only one Eg orbital per site. If we were to include the two degenerate orbitals. 



we would need to consider also the Coulomb and exchange interactions between them to 
produce the Mott insulating states at both ends of the concentration range. S and a are the 
Pauli matrices for spin | at site i for localized and itinerant electrons respectively, is the 
on-site energy, t the hopping parameter between nearest neighbors, U the on-site Coulomb 
repulsion between two itinerant electrons, and J is the ferromagnetic (J > 0) coupling 
between the localized and itinerant electrons. This Hamiltonian represents each Mn^"*" ion 
at site 2 by a spin Si = |, on which one electron can be added to produce Mn^"^ . When 
an electron is added in the d-shell of site i, an exchange coupling J is included to favor 



parallel alignment of the added electron to the already existing spin Also to avoid the 
possibility of Mn^+ we include a strong Coulomb repulsion U and we take U oo. Without 
losing essential physics we simplify further by taking only the z component of the exchange 
interaction. Thus the states of the system are characterized by itinerant electrons moving 
on a frozen distribution of localized up or down spins. To obtain site Green functions and 
thus local density of states for this problem, we ignore at the start the site dependence of 
the diagonal energies: i.e. we set Cj = e and we use an alloy analogy approximation to obtain 
the effect of J (assumed larger than t) in the electronic band structure of the model. Using 
the Renormalized Perturbation Expansion |]T3[ in the manner described in |]rT| we obtain 
the corresponding local Green functions and the average density of states for spin up and 
down. The densities of states for each spin split into two bands centered at E± = (e ± J) 
with weights and widths that depend on the number of sites with each spin Sz = +1/2 or 
—1/2. i.e. they depend on the magnetization of the system. The electronic structure of the 
compounds consists of essentially four bands, two for spin up and two for spin down. The 
splitting between the up and down bands is given by the intra-atomic exchange energy J, 
their weight and width by the normalized magnetization m = 2 < S >. The Fermi level falls 
always in the lower bands so that the transport properties are determined by these bands. 
Consequently, for J » \/Kt , where K is the connectivity, using the site density of state 
(Eq. (11) in Ref. W^) the averaged density of states per site reduces to 
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z/^(ir + l)y'4irtV^-(^-i?)2 

where E = {e — J) and = {1 + fi m)/2 (/x = ± for up and down spin respectively). 

At this point, we introduce the effect of the disorder originated by the substitution of 
some of the rare earth ions by Sr,Ba or Ca. We assume that this can be described within the 
model by making the diagonal energies site dependent. As is well known, since Anderson's 



original paper |]T6| a distribution of diagonal energies produces localization of the electronic 
states from the edges of the bands to an energy within them which is called " mobility edge" 
(ME). The precise position of the ME is difficult to calculate and different localization criteria 



result in different values for it [jT^. However, we do not aim here to an absolute value for 
the ME but rather to its change with respect to the Fermi level when the magnetization 
changes from saturation to cero. For this reason we assume that there is no localization 
before disorder and for simplicity, we use a Lorentzian distribution of energies ( width F 



and the Ziman criterium of localization 19 



From the ensemble-averaged Green function we obtain densities of states. 

/ + 00 
Po,{e')L{e-e')de', (3) 
-oo 

where L{x) is a Lorentz distribution given by 

Within this comparative approach one can make the further approximation of replacing 



in Eq.(3) Pq^ by a square density of states with the same width = 2t^K and the 
same weight to obtain, 

p^{e) = {arctan [{W^ - e)/F] + arctan [{W^ + e)/F]} , (5) 

which allows for analytical expressions for the number of particles n, and the internal 
energy E as functions of the magnetization m , and the Fermi energy ep- In some instances, 
when the Fermi level falls too near the band edge, this approximation can differ from the 
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more realistic case where the density of states increases as ^ . We will see bellow that this 
is the case for n=0.15 in the samples we use to compare our results with. 

To proceed further, we need an expression for the entropy of these system. Again for 
comparative purposes, we resort to the simplest possible form compatible with our earlier 
approximations, that of a spin one half array of sites: 

5 = In (2) - Z/+ In (2z/+) - z/_ In (2z/_) . (6) 

More accurate forms of the entropy valid in the mixed valence regime can be used, see 
for example pO . 



In the presence of a magnetic field if, the free energy per site is then , 

G = E-TS- fiBmH, (7) 

where T is the temperature and /^b is the magnetic moment per site. 
We proceed as follows: for each n, we use (assuming ksT « Wa) 



n 



T.rp,i^)de, (8) 



to obtain a relation between n, m and ep from which ep can be determined numerically. 

The free energy is then a function of m and T only and allows, by minimization, to 
determine m{T). The resulting m{T) (Shown in Fig.l) does not differ essentially from 
the law of corresponding states for spin 1/2. Having obtained m(T) for each value of the 
parameters we can determine the up and down mobility edges ( 5+ and -B_ ) and the Fermi 
Energy. They are also plotted as functions of temperature in Fig.l. 

Following Mott and Davies |]2T| we calculate the transport properties assuming that two 
forms of d.c. conduction are possible: thermally activated hopping and excitation to the 
mobility edge. When the difference between the Fermi level and the mobility edge A is not 
too large as compared to ksT , the conductivity is dominated by particles in the extended 
states, and is given by the usual relaxation time form. 



(9) 



in which a is the Mn-Mn distance in the simple cubic lattice, /(e) is the Fermi function. 
We assume that the relaxation time is a step function equal to zero for e < and takes 
a value Tq related to the minimum metallic conductivity for e> ,where according to Ref. 
ITSl B^ = —Jt'^K'^ z/^ — . Further replace v^{e) by its average = W^ia? /2h^ to obtain: 




An Anderson transition takes place when B^ vanishes. For {t^K"^ — F^) < all 
eigenstates became localized. 

III. RESULTS AND DISCUSSION 

In what follows we take K = 5 appropriate to describe the simple cubic lattice of the Mn 
sites and t = 1 fixes the scale of energies. As a consequence of the structure of the model 
and of the approximations that led us to this point, the model becomes symmetric under 
electron -hole transformation in the lower spin up and down bands. 

For n = 0.5 , the Fermi energy vanishes independently of the value of the magnetization 
and one can obtain an analytical expression for the free energy, from which we derive Tc '■ 

Tc = [(F^ + SOt^) arctan(v^t/F) - vVlOt]/ (SnVlOt). (11) 

Connected to the transport properties we can define a characteristic temperature Tm at 
which the mobility edge crosses the Fermi level. Notice however that this crossing does not 
imply any discontinuous change in the resistivity, the only non-analyticity occurs at Tc. For 
n = 0.5 we obtain an explicit expression for Tm '■ 

= n^i^ ^'"[(^^^ ~ arctan(A,/r) - 2A,rin(4 + F^)], (12) 

where = 2t^jKrn^ , = (1 + /i mc)/2 , and rric = — 1. In Fig. 2 we show Tc 
and Tm as a function of F for n = 0.5. 

In what follows we consider n < 0.5 and identify n with the number of holes, which 
we take to be equal to the concentration of divalent component of the alloy. We define as 
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insulator the state where the Fermi level falls below the ME ( A = (B+ — ei?) > 0) . So that, 
for small F the Fermi level falls above the ME ( A < 0) and only the metallic state appears. 



When F increases, A reduces and, finally A = for a critical value F_ = yO.bK'^t'^ — 
(where rric = and Tq = Tm)- When F is increased from F_ , Tm reduces and finally 



Tm = at a critical value F+ = y K'^t'^ — e\ . Above F+ the system remains insulating at 
all temperatures. Consequently, only for F_ < F < F_|_ the transition between metallic and 
insulating regimes appears. All these facts are depicted in Fig. 2 for ri = 0.5 {ep = ). 
Note the similarity of Tc vs F with Tc versus electron-phonon coupling in Millis et al. [|lOl 
In Fig. 3 we show Tc and Tm as functions of n for some values of F. As a consequence 
of the density of states being modified by disorder, the Curie temperatures decrease with 
F, while the increase with n is just a consequence of the energetics of the bands. Tentative 



fitting of the calculated resistivity with the data on Lai-ajSr^-Os of reference [0 gives a 
value of 1.8 t for F. 

In Fig. 4 we have tried tentatively to fit the logarithm of the resistivity as obtained 
from Eq.(lO) to the measurements of Tokura et al. [|T3|.We have chosen to compare with 
these samples to avoid the complications that arise from strong coupling to the lattice in 
the smaller radius compounds [^. To do that, we fix arbitrarily the value of F at 1.8 t. 
We let t vary from sample to sample to fit Tc. Starting with the curve corresponding to 
X or ?T,=0.175 we choose t=1704 K and change to t=1529 K for n= 0.2, to t=1216 K for 
n=0.3, and to t=1600 K for n=0.15. These values of t correspond to bandwidths that range 
between 1.3 eV to 0.93 eV. We then multiply the values of each calculated resistivity by a 
constant ( in the logarithmic plot corresponds to shifting the curves up and down ) to fit 
approximately the value at the maximum. This last constant corresponds to different values 
of To in Eq.lO. which range in the 10~^^ to 10~^^ sec. These Tq's correspond to the minimum 
conductivity defined in Mott and Davis |21| . We can see that the fitting is better in the more 
"metallic " samples than in the n=0.15 sample where one could expect the contribution of 
localized states to be larger and the model results differ more from experiment. Indeed, as 
pointed out above, the resistivity calculated with the square density of states differs even 



more from experiment that the one shown in the Fig.4, which is calculated with the more 
reahstic density of states of Eq. 2 . 

We conclude from the comparison that the model allows to characterize the resistivity 
behavior of different samples by two parameters, one associated to the degree of disorder 
(r), and the other to the hopping energy t. The values of the hopping energy t can be 
affected by displacement of the oxygen atoms, or by polaronic or other many body effects. 

In Fig. 5 we show the magnetic field effect on the resistivity and compare again with the 
results obtained in |jl3l. Here again, we take T = 1.8, t = 1789 K and select Tq = 0.96* 10~^^ 
to fit the H = 8T curve. We take /ig = 0.964* 10~^° erg/Gauss to fit the rest of the curves. 
Indeed, the fitting of resistivity curves in the absence of magnetic field should be taken with 
care because of the effect of magnetic domains walls. 

Three main interactions should be incorporated in a more complete description of the 
whole family of "colossal magnetoresistance " Mn perovkites. 

1 . Static and dynamic lattice effects can modify not only the values of both parameters, 
t and r. but also the thermodynamics of the transition, leading to first order transitions as 
those found in many of the compounds [^. The connection to the dynamics of the lattice 



has been recently very elegantly demonstrated by Zhao et al. 

2. Coulomb interactions between ions, that in combination with point 1 above could also 
produce charge ordering and lead to the reentrant behavior found in [^ . 

3.Superexchange interactions between the localized spins, that lead to canted states, as 



those found in Electron doped Cai-ajY^-MnOs p5 



The thermopower can also be calculated in a similar way. We will report results for this 
quantity in a forthcoming paper. Measurements of this quantity and resistivity in the same 
crystalline samples would be highly desirable. 

To summarize, we have shown that a very simple estimation of the effect of disorder on 
the double exchange mechanism allows to understand resistivity and magnetoresistivity of 
Sr doped La manganites. The most natural source of disorder is the substitution of rare 
earth by Sr, Ca or Ba, but polaronic or other many body effects may act in a similar way. 
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Figure 1. Zero field magnetization (upper panel) and A = — ep) (lower panel) are 
plotted as a function of the normalized temperature {T/Tc) for T = 1.8, K = 5, t = 1, and 
different values of doping n. 

Figure 2. Phase diagram for n = 0.5 {ep = ). Ferromagnetic Tc (solid line) and metal- 
to-insulator Tm (dashed line) transition temperatures vs F, for H = 0, K = 5, and t = 1. 
Regions labelled as FMM (ferromagnetic metal: m ^ and A < 0), FMI (ferromagnetic 
insulator: m ^ and A > 0), FMM (paramagnetic metal: m = and A < 0), and FMI 
(paramagnetic insulator: m = and A > 0). Dotted line is a guide to the eye. 

Figure 3. Fhase diagram. Ferromagnetic Tc (solid lines) and metal-to-insulator Tm 
(dashed lines) transition temperatures vs doping n, for H = 0, K = 5, t = 1, and different 
values of F: (a) F = 1.8, (b) F = 3, and (c) F = 4.5. Regions labelled as FMM (ferromagnetic 
metal: m ^ and A < 0), FMI (ferromagnetic insulator: m ^ and A > 0), FMM 
(paramagnetic metal: m = and A < 0), and FMI (paramagnetic insulator: m = and 
A > 0). Dotted lines are a guide to the eye. 

Figure 4. Zero field resistivity (solid lines) on a logarithmic scale vs temperature in Lai_„ 
Sr„ Mn O3 taken from Ref. The dashed lines represent the fits with Eq. (10) for t = 1 
, F = 1.8, K = 5, and the corresponding values of doping: n = 0.15 , n = 0.175, n = 0.2, 
and n = 0.3. 

Figure 5. Magnetoresistance in Lai_„ Sr„ Mn O3 {n = 0.175). The solid lines show the 
temperature dependence of resistivity in magnetic fields taken from Ref. [|13|. The dashed 
lines represent the fits with Eq. (10) for t = 1 , T = 1.75, K = 5, and the corresponding 
values of magnetic field: B = 0,B = 3T,B = 8T, and i? = 15 T. 
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